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We present experimental measurements and theoretical calculations of photoionization time de- 
lays from the 3s and 3p shells in Ar in the photon energy range of 32-42 eV. The experimental 
measurements are performed by interferometry using attosecond pulse trains and the infrared laser 
used for their generation. The theoretical approach includes intershell correlation effects between the 
3s and 3p shells within the framework of the random phase approximation with exchange (RPAE). 
The connection between single-photon ionization and the two-color two-photon ionization process 
used in the measurement is established using the recently developed asymptotic approximation for 
the complex transition amplitudes of laser- assisted photoionization. We compare and discuss the 
theoretical and experimental results especially in the region where strong intershell correlations in 
the 3s — kp channel lead to an induced "Cooper" minimum in the 3s ionization cross-section. 



I. INTRODUCTION 

Attosecond pulses created by harmonic generation in 
gases [1, 2] allow us to study fundamental light-matter 
interaction processes in the time domain. When an ultra- 
short light pulse impinges on an atom, a coherent ultra- 
broadband electron wave packet is created. If the fre- 
quency of the pulse is high enough, the electronic wave 
packet escapes by photoionization [3] . As in ultrafast op- 
tics, the group delay of an outgoing electron wave packet 
can be defined by the energy derivative of the phase of 
the complex photoionization matrix element. When pho- 
toionization can be reduced to one non-interacting angu- 
lar channel (L), this phase is the same as the scattering 
phase 77l, which represents the difference between a free 
continuum wave and that propagating out of the effective 
atomic potential for the L-angular channel. In fact, the 
concept of time delay was already introduced by Wigner 
in 1955 to describe s-wave quantum scattering [4]. In 
collision physics, with both ingoing and outgoing waves 
the (Wigner) time delay is twice the derivative of the 
scattering phase. 

In general, photoionization may involve several 
strongly interacting channels. Only in some special cases, 
the Wigner time delay can be conveniently used to char- 
acterize delay in photoemission. One such case might be 
valence shell photoionization of Ne in the 100 eV range 
[5, 6]. In this case, there is no considerable coupling be- 
tween the 2s ep and 2p es or ed channels and ed 
is strongly dominant over es, following Fano's propensity 
rule [7] . The case of valence shell photoionization of Ar in 
the 40 eV range [8] is more interesting. In this case, the 
3s photoionization is radically modified by strong inter- 
shell correlation with 3p [9]. As a result, the 3s pho- 
toionization cross-section goes through a deep "Cooper" 
minimum at approximately 42 eV photon energy [10]. 
Such a feature is a signature of inter-shell correlation and 
cannot be theoretically described using any independent 
electron, e.g. Hartree-Fock (HF) model. 

Recent experiments [5, 8] reported the first measure- 
ments of delays between photoemission from different 
subshells from rare e^as atoms, thus raising considerable 



interest from the scientific community. Different meth- 
ods for the measurements of time delays were proposed, 
depending on whether single attosecond pulses or at- 
tosecond pulse trains were used. The streaking tech- 
nique consists in recording electron spectra following ion- 
ization of an atom by a single attosecond pulse in the 
presence of a relatively intense infrared (IR) pulse, as 
a function of the delay between the two pulses [11, 12]. 
Temporal information is obtained by comparing streak- 
ing traces from different subshells in an atom [5] or from 
the conduction and valence bands in a solid [13]. On 
the other hand, the so-called RABBIT (Reconstruction 
of attosecond bursts by ionization of two-photon transi- 
tions) method consists in recording photoelectron above- 
threshold-ionization (ATI) spectra following ionization 
of an atom by a train of attosecond pulses and a weak 
IR pulse, at different delays between the two fields [14]. 
Temporal information on photoionization is obtained by 
comparing RABBIT traces from different subshells in an 
atom [8]. The name of the technique, which we will use 
throughout, refers to its original use for the measurement 
of the group delay of attosecond pulses in a train [15]. 

Both methods involve absorption or stimulated emis- 
sion of one or several IR photons, and it is important to 
understand the role of these additional transitions for a 
correct interpretation of the measured photoemission de- 
lays. A temporal delay difference of 21 as was measured 
for the photoionization from the 2s and 2p shells in neon 
using single attosecond pulses of 100 eV central energy 
[5]. Interestingly, the electron issued from 2p shell was 
found to be delayed compared to the more bound 2s elec- 
tron. Similarly, delay differences on the order of ~ 100 as 
were measured for the photoionization from the 3s and 
3p shells in argon using attosecond pulse trains with cen- 
tral energy around 35 eV. Again, the 3p electron appears 
to be delayed relative to the 3s-electron, with a difference 
which depends on the excitation energy [8] . 

These experimental results stimulated several theoret- 
ical investigations, ranging from advanced photoioniza- 
tion calculations, including correlations effects [6], time- 
dependent numerical approaches [5, 16-18] to semi- ana- 
lytical developments aiming at understanding the effect 



2 



of the IR field on the measured time delays [19-21]. The 
picture which is emerging from this productive theoreti- 
cal activity is that when the influence of the IR laser field 
is correctly accounted for, such time delay measurements 
may provide very interesting information on temporal as- 
pects of many-electron dynamics. 

The present work reports theoretical and experimental 
investigation of photoionization in the 3s and 3p shells in 
argon in the 32-42 eV photon energy range. Besides pro- 
viding a more extensive description of the experimental 
and theoretical methods in [8], we improve the results in 
three different ways: 

• We performed more precise measurements using a 
stabilized Mach- Zehnder interferometer [22] for the 
RABBIT method. The stabilization allows us to 
take scans during a longer time and thus to extract 
the phase more precisely. Some differences with the 
previous measurements are found and discussed. 

• For the comparison with theory, we determined the 
phases of the single-photon ionization amplitudes 
using the Random Phase Approximation with Ex- 
change (RPAE) method, which includes intershell 
correlation effects [9, 23, 24]. This represents a 
clear improvement to the calculations presented in 
[8], using Hartree-Fock data [25], especially in the 
region above 40 eV where photoionization of Ar 
passes through an interference minimum, owing to 
3s-3p intershell correlation effects. 

• Finally, we improved our calculation of the phase of 
a two-photon ionization process, thus making a bet- 
ter connection between the experimental measure- 
ments and the single photoionization calculated 
phases [20]. 

The paper is organized as follows. Section II presents 
the experimental setup and results. Section III describes 
the phase of one- and two-photon ionization processes 
using perturbation theory in an independent-electron ap- 
proximation. Section IV includes inter-shell correlation 
using the RPAE method. A comparison between theory 
and experiment is presented in Section V. 



II. EXPERIMENTAL METHOD AND RESULTS 
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FIG. 1: (color online) Schematic illustration of our experi- 
mental setup. 



the pump attosecond pulse train, but reflecting the outer 
portion of the IR probe beam. The exact position of 
the recombination mirror with respect to the focal po- 
sition of the pump arm is essential in order to precisely 
match the wavefronts of the probe and XUV (extreme 
ultraviolet) beams. A toroidal mirror (f = 30 cm) focuses 
both beams into the sensitive region of a magnetic bottle 
electron spectrometer (MBES), where a diffusive gas jet 
provides argon as detection gas. The relative timing be- 
tween the ultrashort IR probe pulses and the attosecond 
pulse train can be reproducibly adjusted on a sub-cycle 
time scale due to an active stabilization of the pump- 
probe interferometer length [22]. 
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The experiments were performed with a Tita- 
nium: Sapphire femtosecond laser system delivering 
pulses of 30 fs (FWHM) duration, centered at 800 nm, 
with 1 kHz repetition rate, and a pulse energy of ~ 3 mJ. 
A beam splitter divides the laser output into the probe 
and the pump arm of a Mach- Zehnder interferometer (see 
Fig. 1). The energy of the probe pulses can be adjusted 
by a A/2-plate followed by an ultra-thin polarizer. The 
pump arm is focused by a f = 50 cm focusing mirror into 
a pulsed argon gas cell, synchronized with the laser repe- 
tition rate, in order to generate an attosecond pulse train 
via high-order harmonic generation. An aluminum fllter 
of 200 nm thickness blocks the fundamental radiation and 
subsequently a chromium fllter of the same thickness se- 
lects photon energies of about 10 eV bandwidth in the 
range of harmonics 21 to 27. 

The probe and the pump arm of the interferometer 
are recombined on a curved holey mirror, transmitting 



FIG. 2: Electron spectrum obtained by ionizing Ar with four 
harmonics of orders 21, 23, 25 and 27. The ionization channels 
are shown on the top. 

Figure 2 presents an electron spectrum obtained by 
ionizing Ar atoms with harmonics selected by both Al 
and Cr fllters, with orders ranging from 21 to 27. We 
can clearly identify three ionization channels towards the 
35 V, 3s^/, and 3s^3p^ni {ni = 4p or 3d) continua [26]. 
The corresponding ionization energies are 15.76, 29.2 and 
- 37.2 eV. Note that the settings of the MBES were here 
chosen to optimize the spectral resolution at low energy. 
The large asymmetric proflle obtained at high electron 
energy can be reduced by optimizing the MBES settings 
differently. The spectrum due to 3p-ionization is strongly 
affected by the behavior of the ionization cross-section in 
this region. The relative intensities of the 21^* to the 27*^ 
harmonics are approximately 0.2:0.7:1:1. 
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FIG. 3: (color online) Electron spectrum as a function of 
time delay between the attosecond pulses and the IR laser. 
The signal strength is indicated by colors. The spectrum on 
the right (3p) follows that on the left (3s) with a factor of 6 
reduction in the color code, and a slight overlap in energy. 



Figure 3 shows a typical RABBIT spectrogram, i.e. 
electron spectra as a function of delay between pump and 
probe pulses. The electron yield is indicated as colors. 
Compared to the spectra obtained with the harmonics 
only. Fig. 3 includes electron peaks at sideband frequen- 
cies, including additional absorption or emission of one 
IR photon (see Fig. 4). The intensity of these sidebands 
oscillates with a delay at a frequency equal to 2a;, a; being 
the IR laser photon energy, according to 



S2q{r) = a + /3cos(2a;r - A(j)2q - A02q) 



(1) 



where a and /3 are constant quantities, independent of 
the delay and 2q represents the total number of IR pho- 
tons involved, i.e. an odd number to create harmonic 
2q — 1 OT 2q — 1 plus or minus one IR photon. A^2q de- 
notes the phase difference between two harmonics with 
order 2q-\-l and 2q — l^ while A^2g arises from the differ- 
ence in phase between the amplitudes of the two inter- 
fering quantum paths leading the same final state [Fig. 4 
(a)], ta = A(j)2q/2uj can be interpreted as the group de- 
lay of the attosecond pulses [15]. We define in a similar 
way r^^^ = A02q/2uj arising from the two-photon ioniza- 
tion process. Since the same harmonic comb is used for 
ionization in the 35 and 3p shells, the influence of the 
attosecond group delay can be subtracted and the delay 
difference r (35) - r {3p) can be deduced. The results 
of these measurements are indicated in Table 1 for side- 
bands 22, 24 and 26. We also indicate in the same table, 
previous results from [8]. It is quite difficult in such an 
experiment to estimate the uncertainty of our measure- 
ment. The stability of the interferometer is measured 
to be ~ 50 as. The relative uncertainty in comparing 
the phase offsets of different sideband oscillations is es- 
timated to be of the same magnitude or even slightly 
better. 

Our measurements agree well with those of [8] for side- 
band 24. For sideband 22, the measurements performed 
in [8] could not resolve the sideband peak from electrons 
ionized by harmonic 27 towards the continuum Ss'^Sp'^ni 
(see Fig. 3). A new analysis done by considering only 
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FIG. 4: (color online) Energy level scheme of the processes 
discussed in the present work; (a) RABBIT principle; (b) 
Different channels in 2-photon ionization from the 3s and 3p 
subshells. 



TABLE I: Time delay measurements 



Sideband 


22 


24 


26 


Photon energy (eV) 


34.1 


37.2 


40.3 


r<2)(3s)-r(2)(3p) 
this work (as) 


-80 


-100 


10 


[8] (as) 


-40 (-90) 


-110 


-80 


Tec (3s) - rcc{3p) (as) 


-150 


-70 


-40 


r^^^(3s) -r^^^(3p)(as) 


70 


-30 


50 



the high energy part of the sideband peak leads to the 
number indicated in parenthesis in the table, which is in 
good agreement with the present measurement. There is, 
however, a difference for the delay measured at sideband 
26. We will comment on this difference in Section V. 



III. THEORY OF ONE AND TWO-PHOTON 
IONIZATION 

To interpret the results presented above, we relate the 
one-photon ionization delays to the delays measured in 
the experiment. Using lowest-order perturbation theory, 
the transition matrix elements in one and two-photon 
ionization are 



-iEct{k\z\i) , 



(2) 



• IE 



Atomic units are used throughout. We choose the quan- 
tization axis (z) to be the (common) polarization vector 
of the two fields. The complex amplitudes of the laser 
and harmonic fields are denoted by E^^ and E^^ with 
photon energies uj and Vt^ respectively. The initial state 
is denoted \i) and the final state |^). The energies of 
the initial and intermediate states are denoted and e^,, 
respectively. The sum in M^^^ is performed over all pos- 
sible intermediate states \v) in the discrete and contin- 
uum spectrum. The infinitesimal quantity e is added to 
ensure the correct boundary condition for the ionization 
process, so that the matrix element involves an outgoing 
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photoelectron. The magnitude of the final momentum is 
restricted by energy conservation to e = = 17 + 

for one photon and e = k!^ / 2 = VL -\- oj -\- ei for two-photon 
absorption. The two-photon transition matrix element 
involving emission of a laser photon can be written in 
the same way, with uj replaced by —cj in the energy con- 
servation relation and E^^ replaced by its conjugate. 

The next step consists in separating the angular and 
radial parts of the wavefunctions. The different angular 
channels involved are indicated in Fig. 4(b). We split 
the radial and angular dependence in the initial state as 
{r\i) = Yi.m.{r)Rn-i-{r) and use the partial wave expan- 
sion in the final state 

(r|fc) = (87r)i ^ i''e-'^-^''^Y£M{k)YLM{f)RkL{r). (4) 

L,M 

We perform the spherical integration in Eq. (1) and ob- 
tain 



L = li±l 



-M 



(5) 



where the reduced dipole matrix element is defined as 
T^^\k) = Lk(l I {R,L\r\Rn,k) (6) 



using 3j-symbols and with the notation / = \/2l + 1. The 
reduced matrix element (6) is real. When the dipole tran- 
sition with the increased momentum L = + 1 is domi- 
nant, which is often the case [7], the phase of the complex 
dipole matrix element M^^^ is simply equal to 



arg[M^^)(A:)] = r]L{k) - L7r/2. 



(7) 



(There is also a contribution from the fundamental field 
which we do not write here, as well as trivial phases, e.g 
from the spherical harmonic when M ^ [20] ) . Similarly, 
for two-photon ionization. 



M^^\k)(x e'^^^^^h-^YLuik) (8) 

L = X±l,X = l^±l 
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1 h\ 


I -M 









rrii J 



where 

Tii)(fc)=m(^ lo){ol ) ^^^^H^-^^)- 

(9) 

Here, we have introduced the radial component of the 
perturbed wave function, 

where the sum is performed over the discrete and contin- 
uum spectrum, k denotes the momentum corresponding 
to absorption of one harmonic photon such that the en- 
ergy denominator goes to zero (a^^/2 = + H). The 



summation can be decomposed into three terms, the dis- 
crete sum over states with negative energy, a Cauchy 
principal part integral where the pole has been removed 
(both these terms are real) and a resonant term which 
is purely imaginary. The important conclusion is that 
in contrast to the radial one-photon matrix element, the 
radial two-photon matrix element is a complex quantity. 

To evaluate the phase of this quantity, as explained in 
more details in [20], we approximate Rkhi'f) and p,^\{r) 
by their asymptotic values. We have, for example. 



TTK: 

i 



ln(2^r) . . ttA 
H + r]xW - -TT 

K, Z 



(11) 



This allows us to evaluate analytically the integral 
{RkL\r\p^x) in Eq. (9). We obtain 



arg 



^Vx{i^) - vUk) + 0cc(A:, A^), (12) 



where (pcdk, /^) is the phase associated to a continuum- 
continuum radiative transition resulting from the absorp- 
tion of IR photons in the presence of the Coulomb po- 
tential. It is independent from the characteristics of the 
initial atomic state, in particular its angular momentum. 
An important consequence is that, when inserting the 
asymptotic form Eq. (11) in Eq. (8), the scattering phase 
r]L IS canceled out, so that the total phase will not depend 
on the angular momentum of the final state. In the case 
of a dominant intermediate channel A, the phase of the 
complex two-photon matrix element M^'^\k) is equal to 

arg[M(2)(fc)] =r/A(K) -A7r/2 + </.ec(fc,K). (13) 

It is equal to the one-photon ionization phase towards the 
intermediate state with momentum n and angular mo- 
mentum A plus the additional "continuum-continuum" 
phase. The difference of phase which is measured in the 
experiment is therefore given by 



A(9. 



2q 



Vx{i^>)-Vx{f^<) + (/)cc(^,^>)-</>cc(^,^<). (14) 



where , /^<^ are the momenta corresponding to the 
highest (lowest) continuum state in Fig. 4(a). Dividing 
this formula by 2ct;, we have 



where 



r(2)(fc)=rW(fc)+ree(fc), 

T-W(^k) = '^^(^>) ~^a(^<) 



2a; 



(15) 



(16) 



is a finite difference approximation to the Wigner time 
delay drjx/de and thus reflects the properties of the elec- 
tronic wave packet ionized by one-photon absorption into 
the angular channel A. r^^^ also includes a contribution 
from the IR field which is independent of the angular 
momentum. 



Tcc{k) 



c 

2uj 



(17) 



We refer the reader to [20] for details about how to 
calculate Tcc- Fig. 5 shows Tcc as a function of photon 
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FIG. 5: (color online) Continuum-continuum delay Tcc as a 
function of excitation photon energy for the two subshells 3s 
(red) and 3p (blue) for an IR photon energy of 1.55 eV (800 
nm wavelength). 



energy, for the two subshells 3s and 3p and for the IR 
photon energy uj = 1.55 eV used in the experiment. The 
corresponding difference in delays for the 3s and 3p sub- 
shells is only due to the difference in ionization in en- 
ergy between the two shells (13.5 eV). We also indicate 
in Table I the measurement-induced delays for the three 
considered sidebands. 

The processes discussed in this section can be repre- 
sented graphically by Feynman-Goldstone diagrams dis- 
played Fig. 6(a,b). The straight lines with arrows repre- 
sent electron (arrow pointing up) or hole (arrow pointing 
down) states, respectively. The violet and red wavy lines 
represent interaction with the XUV and IR fields. We 
are neglecting here two-photon processes where the IR 
photon is absorbed first [20]. 
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FIG. 6: (color online) Feynman-Goldstone diagrams repre- 
senting one-photon (a) and two-photon (b) ionization pro- 
cesses, (c) Diagrammatic representation of the RPAE equa- 
tions. The second and third diagrams on the right hand side 
refer to time-forward and time-reversed respectively, (d) Two- 
photon ionization including intershell correlation effects, (e) 
Two-photon ionization with the XUV photon absorbed after 
the IR photon. 



IV. INTERSHELL CORRELATION EFFECTS 

To include inter-shell correlation effects, we use the 
random phase approximation with exchange (RPAE) [9] . 
In this approximation, the dipole matrix element of sin- 
gle photoionization is replaced by a "screened" matrix 
element {k\Z\i)^ which accounts for correlation effects 
between the 35 and 3p subshells. These screened ma- 
trix elements, represented graphically in Fig. 6(c), are 
defined by the self-consistent equation: 



\Z\i) 



|^|z)+^l_im^ 



{u\Z\j){jk\V\m) 



- IS 



{3\z\v)kyAy\3i) 

^2 + e,, - e,- 



(18) 



where i and j are 3s or 3p or vice versa and V = l/ri2 is 
the Coulomb interaction. The sum is performed over the 
discrete as well as continuum spectrum. The Coulomb in- 
teraction matrices {jk\V\iyi) and {uk\V\ij)^ represented 
by dashed lines in Fig. 6(c), describe the so-called time- 
forward and time-reversed correlation processes (Note 
that the time goes upwards in the diagrams). If we re- 
place Z by z in the right term in Eq. (18), we obtain a 
perturbative expansion to the first order in the Coulomb 
interaction. More generally, the use of the self-consistent 
screened matrix elements [Eq.(18)] implies infinite par- 
tial sums over two important classes of so-called "bubble" 
diagrams. Each bubble consists of an electron-hole pair 
uj^ which interacts via l/ri2 with final electron-hole pair 
ki. The energy integration in the time-forward term of 
Eq. (18) (first line) contains a pole and the screened ma- 
trix element acquires an imaginary part and therefore an 
extra phase. For a single dominant channel L, the phase 
of the one-photon matrix element [see Eq.(7)] becomes: 

arg[M(i)(A:)] = r]L{k) + S^k) - (19) 

where ^l(^) = h^kL denotes the additional phase due to 
the correlations accounted within the RPAE. The pho- 
toemission time delay is then determined by the sum of 
two terms: 



.(1) 



dr]L 
de 



d6L 
de 



(20) 



The first term represents the time delay in the indepen- 
dent electron approximation, equal to the derivative of 
the photoelectron scattering phase in the combined field 
of the nucleus and the remaining atomic electrons. The 
second term is the RPAE correction due to inter-shell 
correlation effects. 

We solve the system of integral equations (18) nu- 
merically using the computer code developed by Amu- 
sia and collaborators [27]. The basis of occupied atomic 
states (holes) 35 and 3p is defined by the self-consistent 
HE method [28]. The excited electron states are calcu- 
lated within the frozen-core HE approximation [29]. We 
present some results for one-photon ionization in Fig. 7. 
On the top panel, we show the partial photoionization 
cross-sections from the 3p-state calculated using the HF 
and RPAE approximations (see figure caption). From 
this plot, we see that the 3p kd transition is clearly 
dominant at low photon energies. Intershell correlation 
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FIG. 7: (color online) Top panel: The photoionization cross- 
sections a3p^kd calculated in the HF (blue dotted line) and 
RPAE (red solid line) approximations, are compared with the 
asp^ks cross-sections (HF, green open circles, RPAE, pur- 
ple dashed line). The experimental data for the asp cross- 
section are from Ref. [30]. Middle panel: Photoionization 
cross-sections ass calculated in the HF (blue dotted line) 
and RPAE (red solid line) approximations, are compared with 
experimental data [10]. Bottom panel: Correlation-induced 
phase shifts for the 3s and 3p dipole matrix elements. 



effects are more important for the 3p ks than for the 
3p kd transition. The sum of the two partial cross 
sections calculated with the RPAE correction (red and 
green lines) is very close to the the experimental data 
(solid circles) [30]. The middle panel presents the cal- 
culated cross-section for 3s-ionization and compares it 
to the experimental data from [10]. The RPAE correc- 
tion is here essential to reproduce the behavior of the 
cross-section which, in this spectral region, is a rapid de- 
creasing function of photon energy. 

The bottom panel shows the correlation-induced phase 
shifts S^s^kp cind S^p^kd from the same RPAE calcula- 



tion. We observe that the RPAE phase correction 5sp^kd 
is relatively weak. In contrast, Sss^kp varies signifi- 
cantly with energy, especially near the "Cooper" mini- 
mum. This qualitative difference can be explained by 
a different nature of the correlations in the 3p and 35 
shells. In the 3p case, the correlation takes place mainly 
between the electrons that belong to the same shell with 
not much influence of the inter-shell correlation with 35. 
We confirmed this conclusion by performing a separate 
set of RPAE calculations with only the 3p shell included. 
These calculations lead to essentially the same results for 
3p ionization as the complete calculations. In the case of 
intrashell correlation, the time-forward process [see Fig. 
6(c)] is effectively accounted for by calculating the pho- 
toelectron wave function in the field of a singly charged 
ion. It is therefore excluded from Eq. (18) to avoid double 
count. The remaining time-reversed term [second line in 
Eq. (18)] does not contain any poles and therefore does 
not contribute to an additional phase to the correspond- 
ing dipole matrix element. The small phase Ssp^kd is 
due to intershell correlation which is indeed weak. In 
contrast, 35-ionization is strongly affected by correlation 
with the 3p shell. Consequently, the RPAE phase correc- 
tion Sss^kp^ which comes from the correlation with the 
3p shell in the time-forward process, is large and exhibits 
a rapid variation with energy (a tt phase change) in the 
region where the cross section decreases significantly. 

Finally, we generalize our theoretical derivation of two- 
photon ionization to including the effect of inter-shell 
correlation on the XUV photon absorption. As shown 
graphically in the diagram in Fig. 6 (d), we replace the 
(real) transition matrix element corresponding to one- 
XUV photon absorption by a (complex) screened matrix 
element, with an additional phase term. As a conse- 
quence the phase of the two-photon matrix element be- 
comes: 



(21) 



arg[M(2)(fe)] = r]x{^) + 5x{^) - + ^,,{k, k). 

The time delay measured in the experiment is expressed 
as before as r^'^\k) = r^'^\k) ^Tcc{k), with T^'^\k) mod- 
ified by intershell correlation: 



.(1) 



{k) 



2uj 



2uj 



(22) 



Figure 8 presents calculated time delays r^^^ for 3s 
kp^ 3p ks and 3p kd channels. The ionization de- 
lays from the 3p channel do not vary much with photon 
energy and remain small. The 3p ks delay is negli- 
gible while it takes about 70 as more time for the wave 
packet to escape towards the d channel, due to the an- 
gular momentum barrier. The wave packet emitted from 
the 3s channel takes considerably more time to escape, 
especially in the region above 40 eV, owing to strong 
intershell correlation leading to screening by the 3p elec- 
trons. 



V. COMPARISON BETWEEN THEORY AND 
EXPERIMENT 

We present in Fig. 9 a comparison between our exper- 
imental results (see Table 1) and our calculations. The 
dashed blue and solid red lines refer to the independent- 
electron HF and RPAE calculations, respectively. The 
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FIG. 8: (color online) Ionization delay for the three angular 
channels 3p kd (blue dashed), 3p ks (green dashed) and 
3s kp (red solid line). 

circles refer to the results of [8], while the other sym- 
bols (with error bars both in central energy and delay) 
are the results obtained in the present work. Regard- 
ing the two sets of experimental results, they agree very 
well, except for that obtained at the highest energy cor- 
responding to the sideband 26. Our interpretation is that 
we may be approaching the rapidly varying feature due 
to 35 — 3p intershell correlation, leading to a large dis- 
persion of the possible time delays (Note, however, that 
its position seems to be shifted compared to the RPAE 
result). The experimental and RPAE results agree well 
for the first sideband but less for the two higher energy 
sidebands. Surprisingly and perhaps accidentally the HF 
calculation gives there a closer agreement with the ex- 
periment. 



We now discuss possible reasons for the discrepancy. 
Our calculation of the influence of the dressing by the IR 
laser field is approximate. It only uses the asymptotic 
form of the continuum wave functions (both in the final 
and intermediate states) , thereby neglecting the effect of 
the core. This approximation should be tested against 
theoretical calculations, and especially in a region where 
correlation effects are important. We also neglect the in- 
fluence of the two-photon processes where the IR photon 
is absorbed (or emitted) first [31] [see Fig. 6 (e)]. The 
corresponding matrix elements are usually small, except 
possibly close to a minimum of the cross section, where 
the other process, usually dominant, is strongly reduced. 
Interestingly, in such a scenario, the IR radiation would 
not simply be a probe used for the measurement of the 
phase of a one-photon process, but would modify (con- 
trol) the dynamics of the photoemission on an attosecond 
time scale. Finally, in our theoretical calculation, corre- 
lation effects are only accounted for in the single ioniza- 
tion process (XUV absorption). Additional correlation 
effects surrounding the probing, e.g. after the IR photon 
is absorbed might play a role. 

In conclusion, the results shown above point out to 
the need for explicit time-dependent calculations, which 
would account for many-electron correlation and include 
not only one-photon but also two-photon ionization. We 
also plan to repeat these experimental measurements us- 
ing attosecond pulses with a large and tunable band- 
width. Our results demonstrate the potential of the ex- 
perimental tools using single attosecond pulses [5] or at- 
tosecond pulse trains [8]. These tools now enable one to 
measure atomic and molecular transitions, more specif- 
ically, quantum phases and phase variation, i.e. group 
delays, which could not be measured previously. 
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FIG. 9: (color online) Comparison between our theoretical 
calculations (dashed blue line, HF, red line, RPAE) and ex- 
periments (circles [8], crosses, present work) 
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